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Abstract. We present VLA radio continuum measurements of SN1993J in M81 at the frequencies of 0.32 (P- 
band), 1.3 and 1.7 (L-band), 4.9 (C-band), 8.5 (X-band), and 14.9 (U-band) GHz carried out on December 
17 and 21, 2000, about 2820 days after the supernova explosion. We find that a power-law spectrum, free-free 
absorbed by an homogeneous, or clumpy, distribution of ionized gas yields the best fit to the radio data. A 
combined homogeneous- clumpy model is not favored, but neither totally excluded. This result contrasts with the 
modeling of the early (t < 230 days) radio emission from SN1993J, where a mixture of homogeneous and clumpy 
absorbers appeared to be necessary to adequately describe the behavior of the light curves. The radio spectrum 
of supernova SN1993J between 0.32 and 14.9 GHz is well characterized by a — —0.67 ± 0.02 (S v oc typical 
of an optically thin radio supernova. A fit to the radio spectra of SN1993J from ~ 70 up to 2820 days shows 
that the observed spectral index of SN1993J has been slowly evolving since t ~1000 days on, with the observed 
spectral index changing from a ~ — 1 to a = —0.67. This spectral evolution seems to suggest that, in addition 
to the radiative (synchrotron) losses, adiabatic cooling and ionization (Coulomb) losses at the lowest frequencies 
might be contributing significantly to the integrated electron spectrum. 
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The most widely accepted model for the radio emission 
from supernovae is the standard interaction model, in 
which the observed radio (synchrotron) emission is as- 
sociated with the interaction region between the super- 
nova ejecta and the pre-supernova wind of the progenitor 
star. The radio emission shows an early rise and a late de- 
cline due to a low-frequency absorption process. Chevalier 
( 1982 ) considered a number of mechanisms that could ac- 
count for the absorption: (i) free-free absorption by ma- 
terial in front of the forward shock front, (ii) free- free 
absorption by material within the emitting region, (iii) 
synchrotron self-absorption, and (iv) the Razin-Tsytovich 
effect. At the beginning of the 1980's, the available ra- 
dio light curves could be well explained by synchrotron 
radio emission that was partially free-free absorbed by 
ionized thermal material in the circumstellar wind of the 
progenitor star (e.g. Weiler et al. 1986). The availability 



Weiler et al. 1989, 199C) provided further support for this 
mechanism. However, as new radio supernovae have been 
detected, the situation has become more complex. For ex- 
ample, the radio light curves of SN1986J showed a differ- 
ent shape in the rise to maximum, which was interpreted 



of more complete light curves for these supernovae (e.g. 
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by Weiler et al. (199C) as evidence for internal free- free 
absorption. These authors suggested that a filamentary, 
clumpy-like structure of thermal absorbers could also be 
the source of the slow radio turn-on. Chevalier ( 1998) ) has 
recently argued that synchrotron self-absorption is likely 
to be the dominant absorption mechanism for some radio 
supernovae, in particular for SN1987A, and for Type lb 
and Type Ic supernovae. Finally, observations of SN1993J 
seem to indicate that more than one absor ption mecha- 
nism plays a role (Fransson and Bjor nsson 1998 ; Perez- 
Torres et al. |200l| ; Mioduszewski et al. |200l|) . 

SN1993J is to date the radio supernova whose evo- 
lution has been monitored in greatest detail. The bright- 
ness and proximity of SN1993J in M81 offered an unprece- 
dented occasion for VLBI studies, almost since the super- 
nova explosion on 28 March 1993. For the first time, a 
shell-like structure in a young radio supernova was dis- 
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covered (Marcaide et al. 1995a), in agreement with the 
circumstellar interaction model. Marcaide et al. (1995b) 
showed that the expansion of SN1993J was self-similar, 
and produced the first movie of an expanding supernova. 
Later on, Marcaide et al. (1997) reported on the decelera- 
tion of the SN1993J expansion, a result recently confirmed 
by Bartel et al. ( 2000| ). The modeling of the radio light 
curves of SN1993J has also given deeper insight into the 
supernova phenomenon. It seems now clear that the pro- 
genitor of SN1993J had a substantial mass-loss rate dur- 
ing the late stages of its evolution (M> 1CT 5 M Q yr _1 , 



e.g. Van Dyk et al. 1994). The power-law density pro- 
file of the circumstellar medium around SN1993J seems 



to be different from the standard on e (Lu ndqvist 1994 ; 
Van Dyk et al. 1994 ; Marc aide et al. 1997 ; Immler et al. 
2001 ; Mio duszew ski et al. |200l| ), although Fransson and 
Bjornsson ( 1998] ) have succeeded in modeling the SN1993J 
radio light curves with an s=2 density profile for the CSM. 

Based on numerical modelings of the available radio 
data for SN1993J, Fran sson and Bjornsson ( 1998 ) and 
Perez-Torres et al. ( 2001 ) made predictions about its late- 
time, low frequency (< 1.4 GHz) radio light behaviour. 
However, observational results of the low-frequency radio 
emission from SN1993J were lacking. Thus, we carried out 
VLA observations between ^322 and ~ 15000 MHz aimed 
at checking the above modeling efforts, as well as a way to 
discern the main absorbing processes acting in SN1993J. 

2. Observations and Data Reduction 

We observed SN1993J using 25 antennas of the VLA in 
its most extended configuration (A, with baseline lengths 
from 0.68 up to 36.4 km) at L (18 cm) and P (90 cm) bands 
from UT 13:55 to UT 15:25 on 17 December 2000, and at 
U (14.9 GHz), X(8.5 GHz), C(4.9 GHz), L(1.3-1.7 GHz), 
and P(0.32 GHz) bands from UT 08:40 to UT 11:40 on 
21 December 2000. Each frequency band was split into 
two intermediate frequencies (IFs). Data processing was 
made using the Astronomical Image Processing System 
(AIPS). We first describe briefly the observations at the 
higher frequencies, and then emphasize and describe in 
detail those at lower frequencies, particularly those in the 
P-band, deemed to be a discriminator of the dominant ab- 
sorption mechanism in SN1993J. Table ^ summarizes our 
observations. The flux densities reported were obtained 
by combining the data from both intermediate frequen- 
cies (IFs) at each frequency band, except at L-band, for 
which the frequency separation was large enough that we 
could use each flux density value independently. 

2.1. High-frequency (U-, X-, and C-band) observations 

We used a standard continuum mode at U, X, and C 
bands, thus covering a bandwidth of 50 MHz per IF and 
frequency band. At all three bands, we observed SN1993J 
phase-referenced to 0917+624. The source 0917+624 also 
served as the phase reference for the system calibration. 
At U- and X-band, 3C286 was used as primary flux cali- 



brator (assumed of constant flux density), while at C-band 
3C147 was used instead. At U-band, the integration times 
on SN1993J and 0917+624 were t 93 j=21 min and T 09 i7 = 
7 min, while at X-band these were 11.5 and 6 min, and at 
C-band 12 and 7 min, respectively. 

2.2. Low-frequency (L- and P-band) observations 

Observations were done in spectral-line mode in order to 
detect the likely presence of radio frequency interference 
(RFI) in the data, and also to avoid bandwidth smearing. 
In addition, this observing mode allows a more careful 
editing. 3C286 was used at both bands as primary flux 
calibrator, as well as bandpass calibrator. 

L-band observations. Each of the two 12.5-MHz IFs at 
L-band (centered at 1.34 and 1.67 MHz, respectively) were 
split into 8 spectral channels. Observations of SN1993J 
( T 93J = H-5 min) were phase-referenced to 0917+624 
( T 09i7 = 10-5 min). 0917+624 was also used as secondary 
flux calibrator at L-band. Our final hybrid map at L-band 
(not shown here) was made at an effective observing fre- 
quency of 1.50 GHz. The peak of brightness of the map 
(97.4 mJy beam^ 1 ) corresponds to the nucleus of M81, 
and the flux density of SN1993J is 31.9 mJy. The rms of 
the image background is 0.3 mJy/beam. The flux density 
values reported at L-band in Table |l] correspond to each 
IF, 1.34 and 1.67 GHz, respectively. Since the center fre- 
quencies in L-band are quite separated, they can be used 
as independent values for the purpose of fitting a spec- 
tral index (see below). We note that, though the distance 
between SN1993J and M81 (2' 48") is much larger than 
the dimensions of the restoring beam (1.84 x 0.99 arcsec), 
wide field mapping is necessary if one aims at accurate 
determinations of the flux densities of SN1993J at low fre- 
quencies. 

P-band observations. At P-band, each of the two 3.125- 
MHz IFs (centered at 322 and 327 MHz, respectively) was 
split into 32 spectral channels, which permitted a detailed 
inspection of the data, looking for the presence of RFI. 
Since 0917+624 is a relatively weak source at P-band, 
unlike at higher frequencies, confusion of many strong 
sources in the field made it unsuitable as secondary cali- 
brator. Therefore, the much stronger source J1206+6414 
(3C268.3) was used instead. The total integration times 
were T93J «78 min, and T3C268.3 ~39 min. We imaged 
3C268.3 to obtain an accurate value of its flux density. We 
noticed that the coordinates used for 3C268.3 at the corre- 
lator (a = 12 h 06 ro 24!711, S = +64° 14' 40'.'939; J2000.0) 
were off from the actual ones by more than one arcmin. 
Indeed, the coordinates obtained using the AIPS task 
JMFIT were a = 12 h 06 m 24?8978, S = +64° 13' 36'.'425 
(J2000.0). Such position offset for 3C268.3 translated into 
a similar one for SN 1993J and the whole sky around it. 
We corrected the phases of 3C268.3 for this offset. 

With the 3C268.3 data corrected, we imaged SN1993J 
at P-band, aimed to determine its flux density. Since the 
primary beam of the VLA at P-band is so large (« 2° 20'), 
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Fig. 1. (a) Field of of view of the VLA in A configuration (VLA-A) around SN1993J at P-band, on 17 and 21 
December, 2000. The image shows only the inner 2x2 deg (rather than the 4x4 deg actually mapped), which 
includes 38 fields out of the 78 used in the hybrid mapping process. The image is centered on SN1993J. Only 
contours above 20 mJy beam -1 have been drawn, (b) Hybrid map of the galaxy M 82. The contours are (- 
2,2,3,5,7,10,15,20,30,40,50,70,90,100,110,120) x 4.2 mJy beam" 1 , the off-source rms noise. The peak of brightness 
of the map (a = 09 h 55 m 51*814, 6 = +69° 40' 45'.'20; J2000.0) is 537.8 mJy beam -1 . The flux density in the image 
above the 2-rms level is 10.49 Jy. The dimensions of the restoring beam are 8.6 x 7.2 arcsec, with the beam's major 
axis oriented along P. A. =-30. 2 deg. (c) Hybrid map of SN1993J and the nucleus of its host galaxy M81. The contours 
are (-3,3,5,10,15,20,25,30) x 3.3 mJy beam -1 , the off-source rms noise. The peak of brightness of the map corresponds 
to the supernova and is 71.1 mJybeam -1 . The dimensions of the restoring beam are 8.6 x 7.2 arcsec, with the beam's 
major axis oriented along P.A.=-30.9 deg. 



bright sources far from the supernova could greatly affect 
our image. Further, to obtain a good signal-to-noise ratio 
(SNR) in the image, one is forced to map even distant, 
strong sources. The situation is similar to that encoun- 
tered at L-band, where self-calibration was mandatory to 
improve the SNR of the image, but useless unless M81 was 
also included in the imaging process. At P-band, not only 
the M81 nucleus is important, but also far away, strong 
sources such as M82 need also to be taken into account 
to be able to apply self-calibration to the data. We used 
the AIPS task SETFC to create a list of fields with sus- 
pected sources in them, rather than making an enormous, 
computationally unmanageable image. It is worth notic- 
ing that the interferometric phases are initially so poorly 
calibrated that the reconstructed image will show very 
few sources, making very difficult the convergence of the 



hybrid mapping procedure. Therefore, we used the AIPS 
procedure FACES, which searchs for sources in the NVSS 
catalog, to supply an initial model of the sky, and thus 
improve the global phase calibration. We then applied 
wide- field imaging on this self-calibrated dataset: since 
the array configuration looks different to sources in dif- 
ferent parts of the (large) primary beam, we handled this 
by computing the uv-coverage differently for many small 
fields within the task IMAGR. Sources above a thresh- 
old of 33 mJy (at a frequency of 1.4 GHz) were searched 
for as "a priori" candidates in the NVSS catalog, within 
a 2x2 deg region around SN1993J. This resulted in 78 
fields that were cleaned using wide-field mapping. Fig. 2.2 
shows the inner lxl deg of our final, composite image 
of the VLA-A field at P-band. For clarity, only contours 
above 20 mJybeam -1 have been drawn. Panel (b) shows 
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the galaxy M 82, offset about 0f40 from SN1993J. The 
latter is shown in panel (c) along with the nucleus of its 
host galaxy M81. M82 has a flux of ~10.5 Jy, at an epoch 
when SN1993J has a flux of a mere 71.1 mjy (see Table [j] 
for details of the observations). 



unphysical- value for or K3. Clearly, the combined 
"homogeneous-clumpy" model cannot be preferred over 
models 2 or 3, but we also feel we cannot exclude it solely 
on the basis of a 50% increase of the reduced x 2 . 



3. The radio spectrum of SN1993J 
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We fitted the radio spectrum of SN1993J obtained on 17 
and 21 December 2000 using three models for its radio 
emission. These models are described by the equations 
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tions represent a simple power-law spectrum (model 1), 
a power-law spectrum free-free absorbed by a screen of 
homogeneously distributed ionized gas (model 2), and 
a power-law spectrum free-free absorbed by a "clumpy" 
medium (model 3). To facilitate comparisons with the 
modeling of Van Dyk et al. (1994), we have used here 
their notation: is the frequency, in units of 5 GHz; a 
is the observed spectral index; K\ is the flux density at 
5 GHz, in mjy; and K% and Ka are the optical depths 
at 5 GHz of the homogeneous and "clumpy" external ab- 
sorbing media, respectively, at the epoch of our observa- 
tions. Table summarizes our results, and Fig. ^ shows 
the best fits to models 1 through 3. A simple power-law 
spectrum (model 1) yields a good fit to the data from 
14.94 GHz down to 1.34 GHz, but departs significantly 
from the actual flux density at frequencies below 400 MHz, 
and can be ruled out. The best fits are given by models 
2 (homogeneous absorbing medium) and 3 (clumpy ab- 
sorbing medium), which yield an observed spectral index 
a = —0.67 ±0.02. The frequency at which the free-free op- 
tical depth becomes unity is very similar for both model 
2 (v T=1 = 169 MHz) and model 3 {y r=1 = 241 MHz). 
Note, however, that the very-low radio frequency behav- 
ior of the spectrum is significantly different (see Fig. ^). 
While model 2 (solid line) predicts a sharp fall-off below 
frequencies of —200 MHz, model 3 (dashed line) predicts 
a much milder flux density decrease. Unfortunately, the 
absence of radio measurements below 322 MHz does not 
allow us to rule out either of these models. 

Thus, from Table || we cannot exclude, or prefer, any 
of the two alternative models. 

To model the early-time (t < 230 days) radio light 



curves of SN1993J, Van Dyk et al. (1994) invoked a com- 
bined "homogeneous-clumpy" model of the form S v oc 
v§ e~ T (1 — e~ T )/t'. For such model, we can also obtain 
acceptable solutions (e.g. a = —0.67; K\ = 14.4; K 2 = 
8 x 10~ 4 ; Ks = 10~ 5 ) at the expense of increasing x 2 ed by 
w 50%. If the fitting parameters are left unconstrained, 
this model yields an even better fit than models 2 and 
3, but at the unaffordable price of a negative -hence 
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model 3 




Frequency (GHz) 



Fig. 2. Fits to the continuum radio spectrum of SN1993J 
on 17 and 21 December 2000. As indicated in the figure, 
the dotted, solid, and dashed lines correspond to models 
1, 2, and 3, respectively, as described in the text. Note 
that the size of each data point but the one at 0.324 GHz, 
is larger than its associated error bar. 



For a power-law relativistic electron distribution, 
Ne = No E ' , the spectral index of the electron distribu- 
tion 7 is related to the observed spectral index: 7 = 1 — 2 a. 
From our value of a, it then follows that 7 = 2.34 ± 0.04. 
This value of 7 agre es well with that predicted by Fransson 
& Bjornsson ( 199^ ) for t ;> 1000 days. We should note 
that synchrotron self-absorption could play a role in the 
radio emissi on fro m SN1993J (Chev alier 
& Bjornsson |l99|, Perez- Torres et al. [2001 > 
et al. 200 1|) . In particular, the modeling of Fransson 



1998 , Fransson 
Mioduszcwski 



Bjornsson (see their Fig. 12), predicted a flux density of 
~ 70 mjy at 0.32 GHz and for the approximate epoch of 
our observations, while the modeling of Perez-Torres et al. 
predicted at that frequency a flux density larger than 100 
mjy. Such an excess in flux of the latter modeling is likely 
related to the different spectral evolution obtained for the 
electron spectrum. 

Fig. U shows the evolution of the observed spectral 
index of SN1993J from -70 up to 2820 days after the 
explosion. For each epoch, we fitted the available radio 
continuum data (Van Dyk, private communication) us- 
ing model 2 as described above. While at early epochs 
the fits are not superb, they seem to be compatible with 
a k, —1. However, for epochs ;>1000 days, there seems to 
be a clear trend: the observed spectral index, a, becomes 
progressively less steep. This spectral evolution could be 
explained assuming significant energy losses, mainly as- 
sociated with synchrotron, as well as Coulomb and ex- 



pansion losses. Fransson and Bjornsson (199S) proposed 
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Source Name Flux density (mjy) 





14.94 GHz 


8.46 GHz 


4.86 GHz 


1.67 GHz 


1.34 GHz 


0.324 GHz 


SN1993J 


6.7 ±0.1 


10.4 ± 0.2 


14.7 ±0.4 


30.5 ±0.4 


33.7 ±0.4 


71.1 ±3.4 


3C286 (1) 


3293 ± 16 


4939 ± 5 




13708 ± 20 


15157 ± 21 


25940 ± 120 


0917 + 624 (2) 


1164 ± 1 


1370 ± 1 


1472 ± 1 


1272 ± 2 


1194 ±3 




3C147 (1) 






7940 ± 10 






25940 ± 120 


3C268.3 (2) 












10040 ± 20 



^ primary flux calibrator, ^ phase and secondary flux calibrator. All flux density values reported here were obtained by 
imaging each source, and the quoted uncertainty for each flux value corresponds to the off-source rms noise in each image. 



Table 2. Radio spectrum fits for the VLA observations of SN1993J on 17 and 21 December 2000 



Model 


Oi 


tfi(mJy) 


K 2 x 10 4 i/ c , 2 (MHz) 


K 3 x 10 4 ^,3 (MHz) 


Xrcd 


1 
2 
3 


— 0.65 ± 0.02 
— 0.67 ± 0.02 
-0.67 ±0.02 


14.4 ±0.4 
14.4 ±0.2 
14.4 ±0.2 


8 ± 3 169l3g 


17 ± 7 24ltf 2 


8.9 
3.6 
3.6 



a is the observed spectral index; K\ is the flux density at a frequency of 5 GHz; K2 and K3 are the optical depths at 5 GHz 
of the homogeneous and "clumpy" external absorbing media, respectively; x? e d is the \ 2 P er degree of freedom. Models 1 
through 3 are explained in the text. (v c ,2 and i/ c ,3 are not fitted, and correspond to the frequencies at which the free-free 
optical depth is unity for the homogeneous and "clumpy" medium, respectively.) 



a model that reproduces reasonably well the observed 
spectral evolution of SN1993J. In their model, a constant 
fraction of the shocked thermal electrons -characterized 
by a constant spectral index- are injected and acceler- 
ated. These electrons lose their energy mainly due to syn- 
chrotron losses, thus steepening the integrated electron 
spectrum. Moreover, these authors predict that, for the 
physical parameters they are considering, radiative (syn- 
chrotron) losses start to become less important at t > 
1000 days at frequencies below ^2.3 GHz, whereas energy 
losses of the electrons due to adiabatic and Coulomb cool- 
ing contribute significantly below 2.3 GHz, which seems 
consistent with the observational results. We note that a 
multi-frequency monitoring of SN1993J at radio frequen- 
cies down to ~ 70 MHz, would likely prove to be very 
useful to further constrain the population of electrons re- 
sponsible for the emission, as well as to better understand 
the late-time radio evolution of SN1993J. 

Fig. |] shows the fits of model 2 to all data points (solid 
line), and to all but the P-band data point (dashed line). 
Because any absorption model will obviously fit by turn- 
ing down sharply below the lowest frequency data point, 
the availability of radio continuum measurements at the 
lowest possible frequencies is important, in particular as 
a radio supernova ages. Had the P-band data not been 
available, we could have wrongly concluded that the peak 
of the SN1993J spectrum is at v ~ 0.8 GHz at a level of 
~35 mjy. Having the P-data point available, we conclude 
that the peak is around, or below 0.3 GHz, and at least 



-0.5 



-1.5 



100 1000 
Time since explosion (days) 

Fig. 3. Evolution of the observed spectral index, a, for 
SN1993J, as obtained from fitting the radio continuum 
data at each epoch to model 2 (see text for details). 



at a level of 70 mjy. (We note that a similar fit to all 
data, but using model 3 rather than model 2, would not 
change this result.) Similarly, the best fit to our SN1993J 
radio data using models 2 or 3 as described above, but not 
including the P-band measurement, yields a = —0.73, a 
somewhat steeper value than is obtained if all data points 
are used, a = —0.67, but within 2a of it. 
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Fig. 4. Fits to the continuum radio spectrum of SN1993J 
on 17 and 21 December 2000. Although both fits corre- 
spond to model 2, in one case (solid line) all radio mea- 
surements are used, while in the other case (dashed line) 
all but the P-band data point are used. This results in a 
turnover frequency much higher than actually observed. 
See text for details. 



4. Summary 

We present VLA radio continuum measurements of 
SN1993J in M81 taken on 17 and 21 December, 2000, 
about 2820 days after the supernova explosion, spanning 
the frequency range from 0.32 to 14.9 GHz. We discuss 
in detail the P-band (322-327 MHz) observations, since 
the task of obtaining the flux density at these frequencies 
is much harder than at higher frequencies, due to con- 
fusion from many nearby and far away, strong sources, 
e.g. M82. We also point out the importance of having a 
radio spectral coverage that extends down to P-band as 
the supernova ages, since otherwise the conclusions can be 
wrong. 

We discuss three models for the late radio emission of 
SN1993J: (i) a simple power-law spectrum, (ii) a power- 
law spectrum free-free absorbed by a screen of homoge- 
neously distributed ionized gas, and (iii) a power-law spec- 
trum free-free absorbed by a "clumpy" medium. We find 
that the best fit to the data is yielded by a power-law 
spectrum free- free absorbed by either a homogeneous, or 
a clumpy, distribution of ionized gas. However, a mix- 
ture of homogeneous and clumpy absorbers, while not 
preferred, cannot be totally ruled out. The radio spec- 
trum between 0.32 and 14.94 GHz is well characterized 
by a = —0.67 ± 0.02 (S v °c is a ), typical of an optically 
thin radio supernova, and is significantly different from 
that obtained at epochs <;1000 days between 1.4 and 14.9 
GHz, which indicates an ongoing spectral evolution in the 
radio emission from SN1993J. Since the spectral index 7 
of the electron distribution (Ne — NqE~~<), is related to 
the observed spectral index by 7 = 1 — 2 a, our value of a 
translates into 7 = 2.34±0.04, which agrees well with that 



predicted by Fransson & Bjornsson (1998) for t ^ 1000 



days. We also fit the available radio continuum data of 
SN1993J for the period from ^70 up to 2820 days since the 
explosion, using a power-law spectrum free-free absorbed 
by a screen of homogeneously distributed ionized gas. The 
fit shows that the observed spectral index of SN1993J has 
been slowly evolving since t ^1000 days on, with a in- 
creasing from a value close to -1 to -0.67. This spectral 
evolution seems to suggest that, in addition to the radia- 
tive (synchrotron) losses, adiabatic cooling and ionization 
(Coulomb) losses at the lowest frequencies might be con- 
tributing significantly to the integrated electron spectrum. 
An extension of the low frequency observations further 
down to ^70 MHz should help to fine tune the physical 
parameters of SN1993J, and discern between the homoge- 
neous or clumpy nature of its circumstellar medium. From 
a technical point of view, such observations will be chal- 
lenging. 
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